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Abstract 

We analytically study the photorefractive Gunn effect in n-GaAs subjected 
to two external laser beams which form a moving interference pattern (MIP) 
in the semiconductor. When the intensity of the spatially independent part 
of the MIP, denoted by Io, is small, the system has a periodic domain train 
(PDT), consistent with the results of linear stability analysis. When Iq is 
large, the space-charge field induced by the MIP will compete with the PDT 
and result in complex dynamics, including driven chaos via quasiper iodic 
route. 
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Spatiotemporal behavior of Gunn domains in semiconductors is an interesting and 
important current research problem |3|-|8|]. Gunn diodes are not only good systems for 
studying nonlinear dynamics, but also potentially important for various applications J/J |9| . 
In 1996-1997 Segev et al. || and Subacius et al. || showed, respectively, that optical waves 
can excite multiple Gunn-domain formation with spatial periodicity, called periodic domain 
trains (PDT), in deep- impurity doped GaAs and semi-insulating GaAs; both groups used 
the well known Kroemer's criterion |l| to determine the number of Gunn domains. In 1998, 
Bonilla et al. [f| numerically studied the same dynamical system as Segev et al; however, 
their findings are mostly inconsistent with the predictions of Segev et al. Bonilla et al. || 
found that the number of Gunn domains in the system is not consistent with the prediction 
of Segev et al. and in many cases the system shows complex or chaotic behavior. Very 
recently, Shiau et al. used a linear stability analysis (LSA) to study a shallow-impurity 
doped GaAs with length L subjected to two external laser beams which form a moving 
interference pattern (MIP) of intensity I(x,t) = + m cos(Kx + Qt)} and predicted 
PDT in the Gunn diode 0]. Here Q is the frequency difference of the two optical waves, 
K = 27r/A is the interference wave number, A is the grating period and is much smaller 
than L, m is the modulation depth, and Iq is the average intensity. With consideration 
of time-delay feedback, the diode will emit high-dimensional chaotic microwaves, which 
is potentially useful for secure microwave communications, memory devices, applications 
involving photorefractive effects, etc. 0. 

The purpose of this Communication is to clarify different findings of ||, [[|, and || 
using a nonlinear analytic method to study the system of 0j8| as L — * oo. In this approach, 
the intensity of the spatially independent part of the MIP, Jo, plays an important role in 
determining the simple dynamics (i.e., PDT) or complex dynamics of Gunn domains. When 
Iq is small enough, the system shows simple dynamics consistent with the result of the LSA, 
and Kroemer's criterion is still valid. When Iq is large, the space-charge field induced by 
the MIP will compete with the PDT and result in complex dynamics. This picture is also 
useful for understanding the results of other groups |3|-|5j discussed above. 
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We consider two optical waves incident on a n-GaAs with a biased voltage V. The photon 
energy is just above the bandgap of GaAs, i.e., 1.42 eV, so that electron-hole pairs can be 
generated by optical excitation. The generation-recombination processes include complete 
thermal ionization of the donor, generation of electron-hole pairs by the MIP with rate g, 
and recombination of electron- hole pairs with rate 7. Therefore we can write the dynamical 
equations for the electric field E, the free electron density n, and the free hole density p as: 
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Eq. (HD is for the Gauss law. Eq. @ and Eq. (|3]) represent the continuity equations for the 
electrons and holes, respectively. The circuit condition is given at Eq. (|]). Here v(E) dis- 
plays an N-shaped negative differential mobility (NDM). The homogeneous solution of Eqs. 



(|)-(!) is E = V/L, n Q = (N* D + ^/Nff + 4gI /j)/2, and p = (-N* D + y/Ng + 4^ / 7 )/2. 
When I is very small, n w + g Iq / 'N^j , po ~ glo/N^j, and a natural expansion param- 
eter for high order corrections is e = Po/Njr, = gIo/N^ 2/ y. We can rearrange the model equa- 
tions by the following rescaling process for the dynamical variables and parameters in Eqs. 
(H|)- (|): E — > E, n — > n, p — > ep, I —*■ El. The small parameter e can distinguish the 
large variables (i.e., E and n) from small variable p in our model. Rescaling variables in Eqs. 
(|)-(|), we have § dE/dx = e(n-N£)/e+0(e), dn/dt = -d[nv(E)-D n dn/dx]/dx+0(e). 
These results and Eq. @ implies 
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where Jtotif) is time-dependent current density. The physical meaning of Eq. (||) is quite 
clear, which represents the Gunn-domain formation and with drift velocity v(E Q ) [jXJ] . Then, 
Eq. (|5|) confirms our LSA results in Ref. [[!]]. In the following we want to show the detailed 
structure of PDT. 

In order to get the solution of Eq. (|5|), an ansatz of the electric field, we set 
E(x,t) = E + [E s (x,t)e ik ^- V0 ^ + E*(x,t)e- iR( - x - V0 ^}/2, where v = v(E ). E s (x,t) 
is a new variable corresponding to the amplitude of electric- field domains and E*(x,t) 
is its complex conjugate; K (= 2-ajV) is a bulk property of semiconductor. Since K 
is small, the spatiotemporal behavior of the dynamical system within the linear dimen- 
sion A is dominated by the behavior of E s (x,t). Actually, K is an irrelevent parameter 
which is verified by the numerical works of J3J. When E ^> \E S \, we expand v(E) as: 
v(E) « v + v { o\E-E )+vl?\E-E ) 2 /2 + vl; :i \E-E ) 3 /6, where v { j) = d?v{E) / 'd0\ E=Eo . 
Equation (|5]) and above expansions for E(x,t) and v(E) imply 



—N* D v£ J - D n K )E S + (i2D n K -v )—^- -v^> E 2 s —± + 2\E S 



dt v e u u " J s ' y n uy dx 8 u V dx 1 dx J 

With K -> (as length L -> oo) and £ s ~ ^(t)e iX:c , Eq. (0) implies 

f = -^-igiV>S 3 ) + ^ 2) )|^. (8) 

Here C = (eNjjvjp/e + D n K 2 + iKv ). With #(t) = R(t)expiQ(t), Eq. © implies R(t) = 
-(a/6) 1 /2[ r 2 e -2a i/ ( r 2( 1 _ e -2at) + Q )]i/ 2) = _ lKv<£ ] f R 2 (t')di? , where a = 

eNpVQ'/e + D n K 2 , b = eN^v^/8e, and r is related to initial value of i?(t). When a > 
(i.e., the operating point is within the regime of positive differential mobility PDM), b > 0, 
and t — > oo, i?(t) and 0(t) tend to and —Kv^t, respectively. These results mean that the 
origin on the complex plane of \I/(t) is a stable focus, which also implies electric-field profile in 
material is homogeneous, i.e., E(x,t) = Eq. When a < (i.e., the NDM regime), b > 0, and 
t — > oo, = (—a/6) 1 / 2 and 9(t) = —K(v + clVq /8b)t. This solution describes a circular 
trajectory with counterclockwise angular frequency K(vq + qVq /8b), radius (— a/6) 1//2 , and 



transient response time The transition from stable focus to stable limit cycle is a 

supercritical Hopf bifurcation. Since K — > 0, E s ~ ^{t)e lKx and ^(t) = R(t)expiQ(t) imply 
for a < 0, b > 0, and t — > oo: 

E(x,t) = E + (-a/b) 1/2 cos[K(x - v t - av [ o ] t/8b)], (9) 

which means that the PDT has a spatial period A and a traveling velocity Vq + av^ /8b. 
Since Vq (~ 10 7 cm/s) is much larger than |at>o^/8&|, so we neglect the term av^ /8b below. 



In general, E s may be written as 

5=1 q=l V 
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where a g = eNpv\p / e + D n K 2 q 2 . However, Eq. (JTTJ) needs to be corrected when we consider 
the coupling between spatial modes. This is very important for self-organizing systems. For 
example, consider two spatial modes ^ x e tkx and ^ 2 e t2kx with a 2 > > a\ and b > 0. From 
equations for R(t) and 0(t) , we know that |^i| — ► (— ai/b) 1 ^ 2 and ~~ * when i — > oo. 
Then l^l/l^il wm become zero. This result implies that only unstable spatial modes (i.e., 
a q < 0) have contributions for pattern-forming systems. However, if we combine these two 
modes, i.e., E(x, t) = E + E 2 q =i[$? q (t)e iqKx + ^*{t)e- iqKx }/2, and substitute it into Eq. ©, 
we have 

^l = -(a 1 + iKvo)Vi- E ^^ mi ^ m2 ---^ ms , (11) 

^ = -(a 2 + i2Kv )* 2 - ^mi*™.-*™.. (12) 

where m s = ±1 or ±2, /3 Sjm = eN^v^/e + imsv^^K, and ^* = 0. Based on 
the slaving principle ]nj, we know that the stable mode ^ 2 will be determined by the 
unstable mode ^i, and that ^ 2 , and hence the ratio l^l/l^il = r 2,i{K, Np), w ^ nonzero. 
For simplicity, in the following we assume that r q+ i s (= is independent of q 

and denote the number by r, which can be determined by experiments (see below). From 
E(x, t) = E + [E s (x, t) + E*(x, t))/2 and Eq. ([I 7 ]), for large t we get 



V b J 1 + r l — 2r cos[K (x — v t)\ 

Therefore, the maximum and minimum fields are equal to E max = Eq + (— ai/b) l / 2 {\ — 
r) _1 and E min = E — (— a 1 /6) 1 / 2 (l +r)~ 1 , respectively. According to Pockel's effect, the 
difference in refractive index between E max and E m i n is — n 3 74i(— ai/6) 1 / 2 (l — r 2 ) -1 , where 
n is refractive index of pure material and 741 is the electro-optic coefficient (n = 3.3 and 
7 41 = 1.43 x 10 -10 cm/V for GaAs when the wavelength of the incident light is 1 fim). 
The profile of Eq. fli~3|) , shown in Fig. 1, is strongly r dependent. If r is near 0, E(x,t) 
is approxmately equal to cos[K(x — v t)}. If r is near 1, Eq. (|13|) shows an obvious shock- 
wave structure. In some senses, r is related to the product iV^A in Kroemer's criterion for 
domain formation. To see this, we first consider r — > 1. In this case, (|^ 9 | — l^+iD/l^gl, and 
hence (|a 2 | — |ai|)/|ai|, is much smaller than 1 and there are many excited unstable Fourier 
modes. Thus, —eN^v^ /e will be much larger than D n K 2 and it follows from Eq. (|10|) 
that Kroemer's criterion for sharp domain formation is a\A/vo ~ eN^v^ A/(v e) ^> 1, i.e. 
N^A > — ■u ee~ 1 /' u o 1) (~ 10 12 cur 2 ). On the other hand, in the case r — > 0, —eNpVfr'/e is of 
the same order as D n K 2 and Kroemer's criterion for the appearance of sinusoidal electric- 
field domain is ^A/vq « eN^ 1] A/ (v e) « 1, i.e. N^A « -v Q ee~ l /v^ ] (~ 10 12 cm' 2 ). 



Therefore, Kroemer's criterion appears in Eq. (fy|) via the representation of r parameter. 
To determine the relation between r and iV^A, we can measure (E max — E )/(E Q — E min ) 
(= (1 + r)/(l — r)) at different values of A73A. Then we will find that r is a function of 
NpA. In other words, r can be determined by experiments. 

In order to get hole distribution, we assume p(x,t) = Po + [p s (x,t)e tK ^ x+fJ ' pEot ^ + 
p* s (x, t)e~ iR( - x+l * pEo ^]/2, and substitute such p(x,t) into Eq. to find 

^ = —(D P K 2 + jN* D )p s + D p ^ + (2itfD p + /V^o)^ + S(x, t). (14) 

Here S(x,t) = 2gI e- iR ^ x+ ^ E ^mcon(Kx + fit) + po(//p - je/e){^ + iKE s )e~ iR( - V0+ ^ E °^. 
Equation fli"4]) is a linear partial differential equation with external driving function S(x,t) 
including the MIP and multiple Gunn domains. Thus p s can be considered as a linear 



combination of a homogeneous solution Pg(x, t) and an inhomogeneous solution p l s (x, t). 
Since the characteristic length of p s is equal to A (<C L), we may treat K as zero in Eq. 
(0) and write dp h Jdt = —jN^p^ + D p d 2 p h s / dx 2 + fi p E dp^/dx, 

iNlpi + D p ^ + fi p Eo^ + 2gI mcos(Kx + ttt)+po(v P -^)^- (15) 



dt s dx 2 dx 

It is easy to find solution of p^ and show that p^ 



as t 



e dx 
oo. When E. 



oo r q-li 



ai /b) 1 ' ' 2 exp[iqK(x — v t)] is substituted into Eq. (|T5|), we can find solution 



for p\ . Therefore, the hole distribution p(x, t) = p + Re(pl) is 
gl 2gl m. 



p(x, t) 
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oo qr i-±{ q ( VQ + ^Eq) cos[qK(x - v t)] + {jN^ + q 2 K 2 D p ) sm[qK(x - v t)]} 



D J2 



(16) 



( 1 N* D + q 2 K 2 D p ) 2 + q 2 K 2 (v + fi p E ) 2 

Here A = (D P K 2 + 7 iV£) 2 + (K/j p E - fi) 2 and D = gI K(ij p - 7 e/e)(-a 1 /&) 1 /7 7 iV| ) . The 
pattern profiles of Eqs. (|13D and (0) correspond to simple dynamics. 

When Iq is large and the higher order corrections in Eqs. (|l])-(^P need to be included, 
then we have 



dE/dx - e(n - N* D )/e = -ee/ep {0 \ 



(17) 
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dx 



gI( X ,t)- 1 (N* D + -^-)pW 
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where p^ and E^ are solutions of Eqs. (|I6| ) and (|T3|) , respectively. These results imply 
that up to order 0(e) we have 

1 



d E 1 dP d"^ E 

— + -eN* D v(E) + v(E)— - D n — - ^J„(f) = S(x,t), 



(19) 



where S(x, t) = ee 



p^v(E^) - fipE^pW - (D p + D n )dp( 0) /dx] /e and represents the 
electron- hole coupling. To understand the underlying physics of Eq. ([19]), we follow the 
same procedure as the case of simple dynamics, to get 

1 (e - -\ - rA 



~dt 



^N* D v^ + iKv^>) |^| z ^ + ^ S{x,t)e- tKx dx, 
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(20) 



which is a generalization of Eq. (j8|) and represents a self-sustained oscillator (i.e., PDT) 



under the external driving force eF(t)/e (= J S(x,t)e dx) and 



F(t) 



TV . 71 
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where G = gI /lN* D , G, = 2gI m(D p K 2 + ^N* D )/A, G 2 = -2gI m(K^ p E - to) /A, G 3 = 
-gI K 2 (ti p - 7 e/e)(-a 1 /6) 1 / 2 (t;o + fi P E )/B, G 4 = -gI K(fi p - 7 e/e)(-a 1 /6) 1 / 2 ( 7 iV| ) + 
K 2 D P )/B, and B = jN^N^+K 2 Dp) 2 +K 2 (v +fipEo) 2 . The intrinsic oscillating frequency 
in Eq. (|20|) is Vq/A and the external driving frequencies include Vq/A and Q/2tt. If the 
operating point is in the PDM regime (i.e., at > 0), the last three terms in Eq. (|2lD will 
be zero. This means that there is no Gunn-domain formation in the semiconductor, but 
the external optical waves, i.e., the first two terms on the right-hand side of Eq. (pT|). 
will trigger a space-charge field (SCF) ~ e t{Kx+nt)^ w j 1 j c ] 1 j s we \\ k n0 wn in photorefractive 
material . On the other hand, if the operating point is in the NDM regime, in addition to 
the SCF, the PDT will be generated by the intrinsic electrical instability with the oscillating 
frequency Kvq. The relative strength of the SCF and PDT depends on the strength of 
spatially independent part of the MIP, i.e. Iq. When I Q is small, SCF may be neglected 
and the behavior of the system is dominated by the PDT. When I is large, the SCF will 
compete with the PDT and the system will show complex dynamics. For example, these 
two competing frequencies Q/2tt and Vq/A can lead to quasiperiodic behavior. 

The picture presented above is useful for understanding the results of other groups. 
The system numerically studied by Subacius et al. may be represented by Eqs. ([!])- (|]) 
of this Communication with iV^, = and Q = 0. When the system is illuminated by a 
pulse of light interference field with intensity I(x,t) = I f(t)[l + mcos(i^x)] and m << 1 
(f(t) is the temporal shape of a laser pulse), the temporal part Iof(t) generates background 
electrons in the conduction band and the very small spatially dependent part of the light 
makes a small modulation on the background. The situation is similar to the system of 
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the present manuscript with very small Jo- Fig. 1 in the present Communication is very 
similar as the numerical results of Subacius et ai. On the other hand, Segev et al. || and 
Bonilla et al. |J considered deep-impurity doped GaAs under the externally optical waves, 
i.e., I(x,t) = Io[l + mcos(Kx + Qt)]. In their case the modulation depth m also plays 
a crucial role. For example, Segev et al. considered m as a small parameter. Then the 
dynamical variables can be analytically expanded to Fourier series based on the form of 
I(x, t). Therefore they can get multiple Gunn-domain order by order. Actually the physical 
idea of Segev et al. is the same as Subacius et al.. The spatial period of I(x,t) determines 
the spacing of Gunn domains and Iq (or Iof(t)) generates background electrons. However 
the work of Segev et al. didn't discuss the stability of multiple Gunn-domain, which will lead 
to the travelling domains with the same moving velocity (i.e., Q/K) of I(x, t). Furthermore, 
according to the prediction of Segev et al. the moving velocity of electric-field profile in 
PDM or NDM regime has no any difference. This unsuitable prediction can be addressed 
and verified by Eq. (|2p) . As for the numerical work of Bonilla et ai, they showed complicated 
spatiotemporal behaviors of Gunn domains, e.g., quasiperiodic route to chaos, are mostly 
observed in their simulations. In other words, PDT didn't appear in their system. Therefore 
they concluded that Kroemer's criterion can not give the correct numbers of Gunn domains. 
We think that m = 0.1 used in their study is large enough to generate driven chaos due to the 
competition between SCF and PDT. Therefore both considerations of Kroemer's criterion 
and pattern competition may give correct physics in light-triggered Gunn-domain systems. 

In conclusion, using a nonlinear analytic method we confirm the LSA results of and 
find that the competition between the SCF and the PDT could be a dominant factor to 
determine the spatiotemporal pattern in semiconductor. We hope that the formulation of 
the present Communication can be further verified by numerical calculations or experiments. 

We thank Jonathan Dushoff for a critical reading of the paper. This work was supported 
in part by the National Science Council of the Republic of China (Taiwan) under Contract 
No. NSC 89-2112-M-001-084. 
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Figure Caption 



FIG. 1. The profile of E(x, t) represented by Eq. ( |THj ) for several values of r, where Eq = 10 
kV/cm, (—ax/6) 1 / 2 = 9.8 kV/cm, A = 15 /im, and L = 45 yum. 
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Fig. 1 Shiau and Hu 
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